REMARKS/ARGUMENTS 

The Applicants originally submitted Claims 1-21 in the application. No claims have been 
added, amended or cancelled. Accordingly, Claims 1-21 are pending in the application. 

I. Formal Matters and Objections 

The Examiner has objected to the specification citing several errors and has requested an 
additional figure for understanding of the subject matter of the invention noting that FIGURE 2 is 
not helpful in showing detailed steps of the invention. (Examiner's Action, pages 3-4). In response, 
the Applicants have amended the specification to correct the errors cited by the Examiner. The 
Applicants also point out that in FIGURE 1 the designation C3 represents each block connected to 
the particular block containing C3 (a total of eight blocks) just as CI represents two connected 
blocks and C2 represents two connected blocks. C3, therefore, does represent a more coarse 
description than CI or C2 since C3 represents a larger area. Additionally, the Applicants point out 
that FIGURE 2 (FIGURE 2 A as amended) is not a flow chart to illustrate detailed steps but is a block 
diagram to illustrate parts of a system that perform specific functions and/or procedures. The 
Applicants, nevertheless, have added FIGURE 2B to more fully provide an understanding of the 
invention's subject matter as requested by the Examiner. 

In addition, the Examiner has requested a document that was mentioned in the IDS but was 
not submitted. The requested document is "Introduction to Numerical Analysis," which is a text 
book that was cited by the Applicants as a background reference for numerical analysis. The 
Applicants believe that a copy of this text book does not need to be submitted due to its voluminous 
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nature and since the text book was only cited as a background reference. If the Examiner still 
requires a copy of the text book, the Applicants will gladly supply a copy. 

The Examiner has also asserted that the Applicants are required to amend the disclosure to 
include material incorporated by reference. (Examiner's Action, page 3). As the Examiner is no 
doubt aware, the Applicants are only required to incorporate essential material into the specification. 
The Applicants, however, believe that the material incorporated by reference is not essential and is 
merely background information. More specifically, "GMRES: A Generalized Minimal Residual 
Algorithm for Solving Nonsymmetric Linear Systems" is cited to provide background information 
for Generalized Minimal Residual (GMRES) methods. (Page 16, lines 14-18). In addition, 
"Introduction to Numerical Analysis" and "Preconditioned, Adaptive, Multipole-Accelerated 
Iterative Methods for Three-Dimensional First-Kind Integral Equations of Potential Theory" are 
cited to provide background information concerning numerical analysis and capacitance calculations. 
(Page 20, lines 1-10). Since the cited references provide background information, they do not need 
to be incorporated into the specification. 

The Applicants, therefore, respectfully request that the Examiner remove the objection to the 
specification in view of the above discussion and amendments. 

IL Rejection of Claims 1-21 under 35 U.S.C. §112 

The Examiner has rejected Claims 1-21 under the first and second paragraphs of 35 U.S.C. 
§ 1 1 2 for not being enabled and for being indefinite. Specifically, the Examiner asserts that "integral 
equation formulation" is not enabled or adequately defined in the specification and that "integral 
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equation formulation" is a limitation in all independent claims. (Examiner's Action, pages 5-6). The 
Applicants respectfully disagree. 

"Integral equation formulation" does denote the use for the invention in the preamble of each 
independent claim. The Applicants believe, however, that the specification sufficiently enables 
"integral equation formulation" so that one skilled in the pertinent art can make or use the invention. 
Additionally, the Applicants believe that "integral equation formulation" is adequately defined 
within the specification to particularly point out and distinctly claim the subject matter which the 
Applicants regard as the invention. Moreover, one skilled in the pertinent art, especially in the art 
of IC capacitance calculations, understands the term "integral equation formulation." 

In reference to FIGURE 2, the specification describes the system 200 including the integral 
equation formulator 210, the charge variation function generator 220 and the conductive geometry 
generator 230. The system 200 may employ the components 210, 220, 230, when determining the 
capacitance of an integrated circuit using an iterative linear solution method. The specification 
generally describes the interaction of the components 210, 220, 230, when determining the 
capacitance on pages 11-13. Starting on Page 14, the specification describes in more detail the 
operation of the integral equation formulator 210 in solving the iterative linear solution. The 
specification also describes the operation of the integral equation formulator 210 along with the 
charge variation function generator 220 and the conductive geometry generator 230 in describing the 
iterative linear solution process to provide a more complete picture and show the advantages of the 
invention. (Page 17, line 5 to Page 25, line 17). 

In addition, the specification denotes that the "integral equation formulation" may be a Fast 
Distribution Method (FDM). (Page 4, lines 20-21). Additionally, the specification describes the 
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FDM and compares the FDM to a Fast Multiple Method (FMM). (Page 19, line 6 to Page 20, line 
14). The specification, therefore, sufficiently enables and defines "integral equation formulation" 
through the description and use of the integral equation formulator 210, the FDM and methods of 
analysis well known within the pertinent art. 

Accordingly, the Applicants request that the Examiner remove the rejection of Claims 1 -2 1 
under the first and second paragraphs of 35 U.S. C. § 112 in view of the above discussion. 

III. Rejection of Claims 1 and 8 under 35 U.S.C. §102 

The Examiner has rejected Claims 1 and 8 under 35 U.S.C. §102(e) as being anticipated by 
U.S. Patent 6,397.171 to Belk. In the Examiner's Action, the Examiner asserts that Belk teaches 
each and every element of independent Claims 1 and 8. (Examiner's Action, pages 6-7)^ The 
Applicants respectfully disagree. 

Belk teaches modeling metalization structures by selecting representative sub units and using 
the self and mutual interactions of the sub units as an initial solution to describe all interactions 
between similar metalization sub units in an overall system of metals. (Abstract). The sub units are 
a system of structures including straight polygons, bends and intersections which are decomposed 
from the metalization structure. (Column 6, lines 21). 

Belk does not teach, among other things, a system for generating a representation of charge 
distribution for a given capacitive structure including creating a multidimensional charge variation 
function that is independent of a conductive geometry of the capacitive structure. (Claims 1 and 8). 
As defined in the present specification, a charge variation function is not a representation of a charge 
distribution but is a function used to modify the charge distribution function to reach a desired 
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resolution. (Page 23, lines 4-6). Instead of charge variation functions, Belk merely teaches charge 
distributions that are determined for representative subunits of the metalization structure. (Colunm 
1, lines 55-62). The charge distribution functions on each subunit may be decomposed into 
mathematical functions that capture the differing properties of the components of the charge 
distributions. (Column 12, lines 3 1 -37). Additionally, the charge distributions are dependent on the 
subunits which are decomposed structures from the metalization structure. The charge distribution, 
therefore, is dependent on a conductive geometry of the metalization structure. 

Since Belk does not teach creating a multidimensional charge variation function that is 
independent of a conductive geometry of the capacitive structure, Belk, therefore, does not disclose 
each and every element of the claimed invention associated with independent Claims 1 and 8. 
Accordingly, the Applicants respectfully request the Examiner to withdraw the § 102(e) rejection 
with respect to Claims 1 and 8. 

IV. Rejection of Claims 2-7 and 9-21 under 35 U.S.C. §103 

The Examiner has rejected Claims 2-7 and 19-21 under 35 U.S.C. §103(a) as being 
unpatentable over Belk in view of U.S. Patent 6,175,815 to Statzler, a journal article written by K. 
Nabors (Nabors), U.S. Patent 6,345,235 to Edgecomb, etal (Edgecomb), U.S. Patent 6,351,572 to 
Dufour or a combination of thereof 

As discussed above, Belk does not teach creating a multidimensional charge variation 
function that is independent of a conductive geometry of the capacitive structure as recited in 
independent Claims 1 and 8. Since Claim 15 has analogous claim limitations as in Claims 1 and 8, 
Belk also fails to teach or suggest the Applicants claimed invention as recited in Claim 15. 
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Additionally, Belk does not suggest creating a multidimensional charge variation function that is 
independent of a conductive geometry of the capacitive structure. In contrast, Belk merely teaches 
determining charge distribution by decomposing a metalization structure into subunits. Statzler, 
Nabors, Edgecomb and Dufour also fail to teach or suggest creating a multidimensional charge 
variation function that is independent of a conductive geometry of the capacitive structure. 

As discussed above, Belk fails to teach or suggest all of the elements of the inventions recited 
in independent Claims 1, 8 and 15. Since Statzler, Nabors, Edgecomb and Dufour fail to cure the 
deficiencies of Belk, the Examiner cannot establish a prima facie case of obviousness of dependent 
Claims 2-7 and 9-2 1 , which include the elements of the respective independent claims. Therefore, 
the inventions as stated in Claims 2-7 and 9-21 are not obvious over Belk in view of Statzler, 
Nabors, Edgecomb and Dufour since Belk, Statzler, Nabors, Edgecomb and Dufour, individually or 
in combination with one another, do not teach or suggest all of the claim elements. Accordingly, the 
Applicants respectfully request the Examiner withdraw the 1 03(a) rejection and pass Claims 2-7 and 
9-21 to issue. 
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Conclusion ^ 





In view of the foregoing amendment and remarks, the AppHcants now see all of the Claims 
currently pending in this application to be in condition for allowance and therefore earnestly solicits 
a Notice of Allowance for Claims 1-21. Attached hereto is a marked-up version of the changes made 
to the specification and claims by the current amendment. The attached page is captioned " Version 
with markings to show changes made." 

The Applicants request the Examiner to telephone the undersigned attorney of record at (972) 
480-8800 if such would further or expedite the prosecution of the present, application. 



Respectfully submitted, 



HiTT Gaines & Boisbrun, P.C. 



Mark E. Kelley 
Registration No. 45,857 




Dated: fU^^<^,^^^ > 



P.O. Box 832570 
Richardson, Texas 75083 
(972) 480-8800 
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O. KAPUR 5-10 
VERSION WITH MARKINGS TO SHOW CHANGES MADE 
IN THE SPECIFICATION: 

(1) Kindly replace the paragraph beginning on page 10, line 4 with the following 
paragraph: 

—The geometry of the nets CI, C2, C3 is described hierarchically and captured to a level of 
detail needed to determine an accurate solution. For example, in FIGURE 1, the net CI and the net 
C2 are captured with relatively detailed geometric descriptions. However, the net C3 may be [is] 
captured with relatively coarse geometric descriptions.— 

(2) Kindly replace the paragraph beginning on page 7, line 8 with the following 
paragraphs: 

-FIGURE 2A [2] illustrates a block diagram of a system for determining a capacitance of 
an integrated circuit constructed according to the principles of the present invention; 

FIGURE 2B illustrates a flow chart of a method of determining a charge distribution for a 
net constructed according to the principles of the present invention: — 

(3) Kindly replace the paragraph beginning on page 11, line 5 with the following 
paragraph: 

-Tuming now to FIGURE 2A [2], illustrated is a block diagram of a system for determining 
a capacitance of an IC, generally designated 200, constructed according to the principles of the 
present invention. The system 200 comprises an integral equation formulator 2 1 0, a charge variation 
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function generator 220, and a conductive geometry generator 230. The integral equation formulator 
210 deteraiines the total capacitance of all the nets CI, C2, C3 within the IC.~ 



(4) Kindly replace the paragraphs beginning on page 1 3, line 20 and ending on page 24, 
line 10 with the following paragraphs: 

-In one embodiment of the present invention, the characteristic geometry function for the 
set of net surfaces, Xr^ is a predefined function. In an alternate embodiment of the present invention, 
Xr is calculated from information related to the capacitive structure of the IC. Xr is scalar-valued 
function that has a n-dimensional position as the input. To illustrate the form of Xr, suppose the 
geometry to be described is the unit square in the jcy-plane. The form of Xr for this geometry then 
would be: 



where s(x)=l for 0<x<l and s(x) is zero elsewhere, and d(z) is the standard Dirac delta function. 
Those skilled in the art are familiar with the properties and use of Dirac delta functions. 

The function f, which represents the charge distribution without regard to geometry, is 
evaluated over R, the set of surfaces of the nets CI, C2, C3. The surface integrals over R are 
reformulated as volume integrals via the relation: 




m 



where V is the volume that the net comprises. 
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To solve the iterative linear solution, the integral equation formulator 210 uses an initial 
guess g, for the charge distribution, and a starting subdivision of the geornetry of the net. The 
integral equation formulator 210 also creates the projection matrix Pr from Xr 



The bracketed expression above is the standard inner product for functions: 

{f,g)-[fg a) 

Also in this equation, is a primitive projection matrix where the (i,j) entry oiP^ - j^ijjlk • 
Finally the IjS represent the choice of orthogonal polynomials. In one embodiment of the present 
invention, Legendre polynomials form the system of orthogonal polynomials. One skilled in the 
pertinent art is familiar with and the use of Legendre polynomials. Also, ?\ will represent that part 
of Pr that contains information of the geometry of subdivision b. 

Using the initial guess g for the charge distribution and the starting subdivision of the 
geometry of the net, the integral equation formulator 210 determines charge distributions for a given 
potential on the nets C 1 , C2, C3 . From the charge distributions, the integral equation formulator 2 1 0 
determines the capacitance of the nets CI, C2, C3. 

In the determination of the charge distributions, the iterative linear solution uses the 
following fundamental equation that relates the charge and the potential: 
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where G (r, r') is the Green's function for the geometry of the IC and p(r') is the charge density. 
Those skilled in the art are familiar with the properties and use of Green's functions. 

In one embodiment of the present invention, the fundamental equation relating the charge 
and the potential is expressed as: 



where the \|/ is a potential distribution and f represents the charge distribution without regard to the 
geometry. From this equation, the integral equation formulator 210 can determine the charge 
distribution function f. 

In one embodiment of the present invention, the integral equation formulator 210 uses a 
Krylov method based on a variation of a Generalized Minimal Residual Method (GMRES) for the 
iterative linear solution. Background information concerning GMRES is discussed in GMRES: A 
Generalized Minimal Residual Algorithm for Solving Nonsvmmetric Linear Systems , by Y. Saad 
and M. H. Shultz, SIAM Journal on Scientific and Statistical Computing, 7(3):856-69, 1986 
(incorporated herein by reference). 

GMRES is used to solve for the charge distribution function f in the following equation: 



¥ = MP J 



(6} 



m 



or, equivalently: 



P,MPJ=P,f 



where Mis an operator that converts a charge density to a potential distribution and \|/o is a potential 
which is 1 over an individual net and 0 elsewhere. Those skilled in the pertinent ^ are familiar with 
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the application of GMRES. Nevertheless, below is a description of the process both to give a 
complete picture of the iterative linear solution process and to show the advantages of decoupling 
the charge variation from the geometry. 

The integral equation formulator 21 0 starts iterative linear solution with an initial guess g for 
the charge distribution and an initial geometry. The integral equation formulator 210 also determines 
xV? which represents the geometry of the box b in the subdivision, and P*',^, which is the part of the 
projection matrix Pr which covers the box b (see FIGURE 3A for more information concerning 
geometry subdivisions and boxes). From this information, the integral equation formulator 210 
computes a potential \|/ based upon the initial charge distribution guess g and the current subdivision. 

The integral equation formulator 210 also determines the desired potential vi/q: 

p^f, = Pso m 

where o is the constant function 1 over an individual net and Pg is defined as Pr [is] over an 
individual net S. The integral equation formulator 210 then determines the difference between the 
desired potential \|/o and the potential v|/. This difference, v|/o-Y, becomes first charge variation 
function f,. Note that fj is not a representation of the charge distribution described by f. f, is a 
function used to modify the charge distribution function f to reach a desired resolution. 

The integral equation formulator 210 determines the ratio y/p, where y = llf,l[R and p = lk|/ol[R. 
If the ratio is sufficiently small, the iterative linear solution is complete. 

If the ratio is unacceptably large, the integral equation formulator 210 normalizes fj and 
proceeds with the linear iterative solution. The integral equation formulator 210 then invokes the 
charge variation function generator 220 to create a charge variation function which refines the 
description of the charge distribution. At the beginning of the iterative process the charge variation 
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function generator 220 creates a second charge variation function fj. In general, the formula for the 
creation of the j+1^*^ charge variation function from the j^*^ charge variation function is: 

Where M is an operator that converts a charge density to a potential distribution. The (i j) entry of 
the matrix representation of M is: 

Where l^j is the i* moment in box b. Note that the operation of M involves all boxes, not just box 
b. G(r,r') is the Green's function for the set of net surfaces R. Those skilled in the art are familiar 
with the properties and use of Green's functions. 

The charge variation function generator performs the operation M using a variant of the Fast 
Multipole Method (FMM) algorithm called the Fast Distributed Method (FDM). The FDM differs 
from the FMM in several areas. First, the FDM omits the FMM's initial step of computing multipole 
representations from point charges. FDM omits this step because the input is already in terms of a 
charge distribution. Second, the FDM omits FMM's direct point-to-point interactions. Instead, the 
FDM uses a source charge distribution in a cube. The potential distributions are calculated in all 
neighboring cubes, including the source cube itself. 

Third, the FDM omits FMM's final step of interpolating point potentials from the local 
expansions. The FDM omits the last step, because the desired output is a potential distribution. 
Finally, FDM uses Legendre polynomial expansions for both charge distributions and local 
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expansions. The FMM uses multipole expansions for representing charge distributions and local 
expansions for representing potential distributions. 

The FMM algorithm, Legendre polynomial expansions and their use are well known to those 
skilled in the pertinent art. Background information concerning numerical analysis and capacitance 
calculations is discussed in Introduction to Numerical Analysis , by J. Soer and R. Bulirsh, Springer- 
Verlag 1979 and in Preconditioned, Adaptive, Multipole-Accelerated Iterative Methods for Three- 
Dimensional First-Kind Integral Equations of Potential Theory , by K. Nabors, et al , SIAM Journal 
on Scientific and Statistical Computing, 15(3):713-735, May 1994 (both incorporated herein by 
reference). 

The use of the FDM results in a significant savings in time over both the analytical 
calculation using the explicit form of M and the FMM. The savings are due to the omission of 
several FMM steps and the omission of FMM' s direct point-to-point interactions. 

Once the charge variation function generator 220 has created the charge variation function 
fj+i, the integral equation formulator 210 uses this new function to create the (i, j) entry of a 
Hessenberg matrix: 

for all i<j. 

The integral equation formulator 210 uses the Hessenberg matrix to solve for the new coefficients 
of the charge variation functions to be used to further refine the charge distribution function f 
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After the entries in the Hessenberg matrix are computed, the integral equation formulator 2 1 0 
then orthogonalizes fj+, from the other fjS. The charge variation function generator 220 then creates 
the J) entry of the Hessenberg matrix by the expression: 
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This expression is the standard norm of functions defined by: 



After the creation of this last entry of the Hessenberg matrix, the integral equation formulator 210 
then normalizes fj+,. The orthogonalization and the normalization of fj+, create orthonormal charge 
variation ftmctions. As a result, the integral equation formulator 210 is able to obtain the correct 
charge distribution function f very quickly. 

The integral equation formulator 210 then solves for x, the coefficient vector for the charge 
variations, using least-squares on the equation: 

Hx = 05} 

where H is the Hessenberg matrix generated by the integral equation formulator 210 as explained 
above, e, is the unit vector with the form (1,0,0,0,0, . . .)(with j number of zeros). The vector x 
is the coefficient vector for the charge distributions and y is the norm of the first charge variation 
fiinction. y is defined as: 

/ = I|/.IL- 06} 
-28- 



Since the charge variation and the geometry are decoupled, the representation of the fj will not 
become very large as compared to the previous methods. Using the residual of the least squares, the 
integral equation formulator 210 computes the factor r as the 2-norm of the least squares residual. 

The integral equation formulator 2 1 0 then generates a new guess for the charge distribution 
function f using the equation: 

f = g + t^J< 07} 

The integral equation formulator 210 then determines the ratio r/p. Where r is the 2-norm 
of the least squares residual r and P is the norm of the desired potential (defined earlier as P = ll\|/ol[R). 
If the ratio is within acceptable limits, the charge distribution has converged and the integral equation 
formulator 2 1 0 terminates the iterative linear solution. 

If the ratio is not within acceptable limits, the integral equation formulator 2 1 0 uses the new 
guess for the charge distribution fimction f, calculated above, as the starting point for the next 
iteration of the iterative linear solution. However, before the integral equation formulator 210 
performs the next iteration, the integral equation formulator 210 must determine if the current 
subdivision requires a more refined approximation. 

The integral equation formulator 210 determines if the current subdivision b requires 
refinement by calculating a charge-geometry error e. In one embodiment of the present invention, 
the integral equation formulator calculates the charge-geometry error e'* associated with box b in the 
subdivision using the following equation: 
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i=k 
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Where 1^ is the i^*" moment basis function, f represents the distribution of the charge without any 
reference to the geometry of the IC and fj^ is the component of f along the i**" moment for the charge 
function contained in the subdivision b. The summation of ij between k and N-1 represent the high 
order polynomials. Finally, N is the number of moments used in the expansion of f. 

The integral equation formulator 210 then determines if the charge-geometry error is within 
acceptable limits using a charge-geometry error criterion. The charge-geometry error criterion is 
defined as: 



Where S is the minimal subdivision containing the selected net. 8 is a small constant that determines 
the allowable tolerance of the criterion. Also, \\q^\\^ is the magnitude of the charge-geometry error 
in subdivision b. Note that the error threshold is relative to the charge distribution f only on the 
subdivision of the selected net.~ 

(5) Kindly insert the following paragraphs after the paragraph that ends on 24, line 10: 
- Turning now to FIGURE 2B with continued reference to FIGURE 2 A, illustrated is a flow 
diagram of an embodiment of a method of determining a charge distribution for a net, generally 
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designated 250, constructed according to the principles of the present invention. The method starts 
in a step 255 with an intent to determine a charge distribution. 

An initial charge distribution and geometry are provided in a step 260. The initial charge 
distribution and the initial geometry may be guesses and are used to start an iterative linear solution. 
The initial guess for the charge distribution may be designated g and the initial guess for the 
geometry may be a subdivision of the geometry of the net. 

After an initial charge distribution and geometry are provided, a first charge variation 
fiinction is then determined in a step 265. The first charge variation fiinction may be the difference 
between \|/o and \|/. In one embodiment, the first charge variation fiinction may be determined by 
solving for and v|/ using Equations 6 and 9. 

After determining the first charge variation fiinction, a determination is made if the charge 
distribution function is within an acceptable limit in a first decisional step 270. In one embodiment, 
the acceptance of the charge distribution function may be within an acceptable limit if the ratio y/p 
is sufficiently small. If the charge distribution fiinction is within an acceptable limit, the method 250 
ends in a step 295. 

If the charge distribution function is not within an acceptable limit, the charge variation 
function is refined through an iterative linear method in a step 280. In one embodiment, f | may be 
normalized before proceeding with the linear iterative method. The method 250 creates a charge 
variation fimction which refines the description of the charge distribution employing Equations 1 0- 
IZ 

After refining the charge variation fiinction, a determination is made if the charge distribution 
function is within an acceptable limit in a step 282. In one embodiment, the charge distribution 
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function is within an acceptable limit based on the ration r/p. If the ratio is within acceptable limits, 
the charge distribution has converged and the method 250 continues to step 295 and ends. 

If the charge distribution function is not within acceptable limits, a determination is made 
if the geometry needs refinement in a third decisional step 287. In one embodiment, the method 250 
may employ Equations 18 and 19 to determine if the geometry needs refinement. If the geometry 
does not need refinement, the method continues to step 280. If the geometry does need refinement, 
the geometry is subdivided into subdivisions in a step 290. After subdividing the geometry, the 
method 250 continues to the step 280. 

While the methods disclosed herein have been described and shown with reference to 
particular steps performed in a particular order, it will be understood that these steps may be 
combined, subdivided or reordered to form an equivalent method without departing from the 
teachings of the present invention. Accordingly, unless specifically indicated herein, the order and/or 
the grouping of the steps are not limitations of the present invention. — 
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